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Abstract. The blazar Markarian 501 (Mrk 501) was observed above 100 GeV with the MAGIC Telescope during May, June 
and July 2005. The high sensitivity of the instrument made possible the detection of the source with high significance in each 
of the observing nights. During this observational campaign, the emitted gamma-ray flux from Mkn 501 was found to vary by 
one order of magnitude, and showed a high correlation with spectral changes. Intra-night flux variability was also observed, 
with flux-doubling times of ~ 2 minutes. The data showed a clear evidence of a spectral peak (in the vFv representation) 
during the nights when the gamma-ray activity was highest. The location of this spectral feature was found to be correlated 
with the emitted gamma-ray flux. In these proceedings we discuss some of the results of this unprecedented spectral and 
i-jh ■ temporal analysis of Mrk 501 observations in the very high energy range. 
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OBSERVATION AND RESULTS 

The observations of Mrk 501 in the Very High Energy (VHE) domain were carried out with the Major Atmospheric 
Gamma-ray Imaging Cherenkov (MAGIC) Telescope (see (2J, 13, La |j] for details of MAGIC). The net observation 
time on the source is 31.6 hours (24 nights) between May and July 2005. Quasi-simultaneously to the MAGIC data, 
Mrk 501 was regularly observed with KVA as a part of Tuorla Observatory's blazar monitoring prograrrQ. And we 
also used data taken with the RXTE satellite's All-Sky-Monitor (RXTE/ASM^ 

The overall Light Curve (LC) of Mrk 501 during the MAGIC observation campaign is shown in Fig. Q] The 
observed flux is shown in three energy bands: VHE (0.15GeV-10TeV), X-rays (2keV-10keV), and optical (1.5eV- 
2.5eV) as measured by MAGIC, RXTE/ASM and KVA, respectively. The X-ray and optical fluxes are computed as 
weighted averages using RXTE/ASM and KVA measurements taken simultaneously with the MAGIC observations 
plus/minus a time tolerance of 0.2 days. A smaller time tolerance substantially decreases the number of X-ray points 
that can be used. The 7-ray flux level of the Crab Nebula (lilac-dashed horizontal line in the top plot) is also shown 
in the upper plot of Fig. [T]for comparison. The Crab nebula flux was obtained by applying the very same analysis 
used for Mrk 501 to the MAGIC Crab nebula data taken during December 2005 under observing conditions similar to 
those for Mrk 501. The estimated Crab Nebula flux level is therefore roughly affected by the same systematics as the 
fluxes obtained for Mrk 501. We found F Crab (>0.15TeV)=(3.2±0.1)x 10~ 10 cirrV 1 , thereafter referred to as Crab 
Unit (c.u.). Fig. Q] shows that the VHE flux from Mrk 501 was about 0.5 c.u. during most of the observation nights; 
yet there are significant deviations from this mean 7-ray emission, with large flux variations occurring in consecutive 
nights. 

During the two nights with the highest VHE activity (>3c.u.), namely June 30 and July 9, Mrk 501 clearly showed 
intra-night flux variations. The corresponding LC in the 0.15-10 TeV band is shown in Fig.|2]with a time binning of 
~2 minutes. A constant line fit to the whole LC gives a % 2 /NDF = 47.9/30 (probability p = 2.0 x 10~ 2 ) for the night 
of June 30, and a % 2 /NDF = 80.6/21 (p = 6.4 • 10~ 9 ) for the night of July 9. Therefore, the emission above 150 GeV 
during the two nights is statistically inconsistent with being constant. The burst's amplitude and duration, as well as 
its rise/fall times, were quantified according to the following function 18J]: 
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1 |http : / /users .utu . f i/kani/lm/| . 

2 The data are publicly available at http : //heasarc.gsfc.nasa. gov/xte_weather /| 



F(t)=a 



(1) 



This model parametrizes a flux variation (flare) superposed on a stable emission: F(t) asymptotically tends to a 
when t — > ±°o. The parameter a is the assumed constant flux at the time of the flare (cf. the horizontal black dashed 
lines in Fig. [2j»; ?o is set to the time corresponding to the point with the highest value in the LC; and b,c,d are left free 
to vary. The latter two parameters denote the flux-doubling rise and fall times, respectively. The resulting fit parameters 
are reported in the insets of Fig. [2] showing flux-doubling times of the order of 2 minutes. This is the fastest flux 
variability ever observed from Mrk 501. These bursts were also studied in four non-overlapping energy ranges; 0.15- 
0.25 TeV, 0.25-0.6 TeV, 0.6-1.2 TeV and >1.2 TeV. The main outcome of this study is a) the burst from June 30 is 
significant only in the energy range 0.25-0.6 TeV, while the burst from July 9 is significantly observed in all energy 
ranges; b) in the July 9 burst, the highest energies are delayed by 4 ± 1 minutes with respect to the lowest energies; 
and c) in the July 9 burst, the relative amplitude of the flux variation increases with energy. Further details from these 
studies can be found elsewhere QJJ] _ 

Mrk 501 showed energy-dependent flux variations throughout the entire MAGIC observational campaign. We 
followed the prescription given in [9] to quantify the flux variability by means of the fractional variability parameter 
F var , as a function of energy. The outcome of this calculation is shown in the left-hand plot of Fig. [3} the flux variability 
increases with energy. The same tendency was reported at X-ray frequencies, but with lower values of F var |5||. 

During our observations, the VHE emission of Mrk 501 was very dynamic, showing significant spectral variability 
on a timescale of days. Nevertheless, most of the data are well described by a simple Power-law (PL) function. This 
does not hold for the two flaring nights of June 30 and July 9 which show clear spectral curvature and thus they are fit 
with a log-parabolic (LP) function. 
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Here, Kq is a normalization factor, a is the spectral index at 0.3 TeV, and b is a curvature parameter. 

The right-hand plot of Fig. [3] shows the energy spectra (EBL corrected) of Mrk 501 for the flaring nights (June 
30 and July 9) for three different flux levels, low, medium and high, which contain 12, 8 and 2 ("non-flaring") nights 
respectively. See caption of Fig.[3]for definition of flux levels. The data shows very clearly a hardening of the spectra 
with increasing flux. Besides, the data from the flaring nights suggest the presence of a spectral peak. The location 
of this peak, determined using the fit parameters and errors from the log-parabolic fit, is 0.85 ±0.13 TeV and 0.44 
± 0.08 TeV for June 30 and July 9, respectively. If this peak exists for the data sets low, medium and high, then its 
location is certainly below 0.1 TeV. Therefore, the location of such spectral peak is very probably correlated with the 
source luminosity. 



CONCLUSIONS 

We have undertaken a systematic study of the temporal and spectral variability of the nearby blazar Mrk 501 with 
the MAGIC telescope at energies > 0.1 TeV. During 24 observing nights between May and July 2005, all of which 
yielded significant detections, we measured fluxes and spectra at levels of baseline activity ranging from <0.5 to >1 
c.u.. During two nights, on June 30 and July 9, Mrk 501 underwent into a very active state with a y-ray emission >3 
c.u., and flux-doubling times of ^2 minutes. An overall trend of harder spectra for higher flux was clearly seen. The 
VHE y-ray variability was found to increase with energy, and it is significantly higher than the variability at X-ray 
frequencies. A spectral peak, at a location dependent on source luminosity, was clearly observed during the active 
states. All these features are naturally expected in synchro-self-Compton (SSC) models of blazar VHE emission. See 
(2] for further details from these observations. 
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FIGURE 1. Multi-frequency LC during the MAGIC observations of Mrk 501 (May- July 2005). Top) MAGIC flux above 0.15 
TeV. The Crab flux is also shown for comparison (lilac dashed horizontal line). Middle) RXTE/ASM 2-10 keV flux. Bottom) KVA 
~ 1 .5-2.5 eV flux. Error bars denote lo~ statistical uncertainties. The X-ray/optical data were selected to match the MAGIC data 
within a time window of 0.2 days. 
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FIGURE 2. Integrated-flux LCs of Mrk 501 for the flare nights of June 30 and July 9. Horizontal bars represent the 2-minute 
time bins, and vertical bars denote lo" statistical uncertainties. For comparison, the Crab emission is also shown as a lilac dashed 
horizontal line. The vertical dot-dashed line divides the data into 'stable' (i.e., pre-burst) and 'variable' (i.e., in-burst) emission 
emission. The horizontal black dashed line represents the average of the 'stable' emission. The solid black curve represents the 
best-fit flare model (see eq. [TJ. The insets report the fit parameters and goodness of the fit. The bottom plots show the mean 
background rate during each of the 2-minute bins of the LCs. The insets report the mean background rate during the entire night, 
resulting from a constant fit to the data points. The goodness of such fit is also given. The background rates are constant along the 
entire night. Consequently, the variations seen in the upper panels of the middle and right-hand plot correspond to actual variations 
of the VHE y-ray flux from Mrk501, thus ruling out detector instabilities and/or atmospheric changes. 
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FIGURE 3. Left-hand: Fractional variability parameter, calculated as described in [9]. Vertical bars denote 1(7 uncertainties, 
horizontal bars indicate the width of each energy bin. Right-hand: Energy spectra of Mrk 501 for the flaring nights of June 
30 (black squares) and July 9 (red up-triangles), and for three data sets which group nights according to whether their integral 
flux above 150 GeV, Fi50G e v (measured in Crab Units (c.u.)), was high (1.0 c.u.<Fi50G e v; green down triangles), medium 
(0.5 c.u.<FigoGeV< 1-0 c.u.; blue open circles), or low (Fi50G e v<0-5 c.u.; pink open squares). Vertical bars denote 1<T uncertainties, 
horizontal bars denote energy bins. Lines show best fits using log-parabolic (for flare nights) and power-law (for high/medium/low 
flux levels) functions. The insets report the spectral indices derived from the fit, as well as the goodness of such fit. The spectra are 
corrected for EBL extinction using [6]'s 'Low' EBL model. 



